Ultra-Deep Mid-Infrared Spectroscopy of Luminous Infrared 

Galaxies at z ~ 1 and z ~ 2 



7 



Dario Fadda 1 , Lin Yan 1 , Guilaine Lagache 3 , Anna Sajina 2 , Dieter Lutz 4 , Stijn 

Wuyts 5 , David T. Frayer 1 , Delphine Marcillac 3 , Emeric Le Floc'h 6 , Karina Caputi 7 , Henrik 

W. W. Spoon 8 , Syl vain Veilleux 9 , Andrew Blain 10 , George Helou 1 
1 IP AC, Cali) p ornia Institute of Technology, Pasadena, CA 91125 

2 Haverford College, Haverford, PA 19041 

3 Institut d'Astrophysique Spatiale, F-9 14-05 Orsay, France; Univ. Paris-Sud 11 and CNRS (UMR 8617) 

4 Max- Planck- Institut fur Extraterrestrische Physik, 85741 Garching, Germany 

5 W. M. Keck Postdoctoral Fellow, Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, 

Cambridge, MA 02138 

6 CEA Saclay, Service d Astrophysique, 91191 Gif-sur-Yvette, France 

SUPA Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh 

EH 9 3HJ, UK 

8 Dept. of Astronomy, Cornell University, 219 Space Sciences Building, Ithaca, NY 14853 
9 Dept. of Astronomy, University of Maryland, College Park, MD 20742 
10 Astronomy Department, California Institute of Technology, Pasadena, CA 91125 
f addaOipac . caltech . edu , lyan@ipac . caltech . edu 



ABSTRACT 

We present ultra-deep mid-infrared spectra of 48 infrared-luminous galaxies in the GOODS- 
South field obtained with the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. These 
galaxies are selected among faint infrared sources (0.14-0.5 mJy at 24 /im) in two redshift bins 
(0.76-1.05 and 1.75-2.4) to sample the major contributors to the cosmic infrared background 
at the most active epochs. We estimate redshifts for 92% of the sample using polycyclic aromatic 
(PAH) and Si absorption features obtaining, in particular, 8 new redshifts difficult to measure 
from ground-based observations. Only few of these galaxies (5% at z ~ 1 and 12% at z ~ 2) 
have their total infrared luminosity dominated by emission from active galactic nuclei (AGN). 
The averaged mid-IR spectrum of the z~l luminous infrared galaxies (LIRGs) is a very good 
match to the averaged spectrum of local starbursts. The averaged spectrum of the z ~ 2 ultra- 
luminous infrared galaxies (ULIRGs), because of a deeper Si absorption, is better fitted by the 
averaged spectrum of HH-like local ULIRGs. Combining this sample with other published data, 
we find that 6.2 /im PAH equivalent widths reach a plateau of ~ 1 /im for £24/1™ ^ 10 n L©. At 
higher luminosities, EW^^^m anti-correlates with Z/24/xm- Intriguingly, high-z ULIRGs and sub- 
millimeter galaxies (SMG) lie above the local EW^^ra-E^^jn relationship suggesting that, at 
a given luminosity, high-z ULIRGs have AGN contributions to their dust emission lower than 
those of local counterparts. A quantitative analysis of their morphology shows that most of the 
luminous IR galaxies have morphologies similar to those of IR-quiet galaxies at the same redshift. 

All z ~ 2 ULIRGs of our sample are IR-excess BzK galaxies and most of them have Lfjr/L 

^ ^ 1 — 16QQA 



ratios higher than those of starburst galaxies at a given UV slope. The "IR excess" (jDaddi et al 



20071 ) is mostly due to strong 7.7 /im PAH emission and under-estimation of UV dust extinction. 
On the basis of the AGN-powered continuum measured directly from the mid-IR spectra, 

we estimate an average intrinsic X-ray AGN lu minosity of L^-ioke V = (0.1=b 0.6) x 10 43 erg/s, a 



value substantially lower than the prediction bv lDaddi et al. (J2] 
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Subject headings: infrared: galaxies - galaxies: star formation - galaxies: high-redshifts - galaxies: 
evolution 



1. Introduction 



The Cosmic i nfrared background (CIB) de- 
tecte d by CO BE (jPuget et al.lll996l : iFixsen et al 



19981 ) peaks around 200 /im with energy com- 



parable to the optical/UV background, imply- 
ing that ~50% of the total optical/UV emission 
is absorbed by dust and re-radiated at mid to 
far-in frared. Using stacking methods, Dole et al. 



(2006) quantified that the CIB emission at 70 
and 160 /im is primarily produced by 24 /im 
sources with S24 Mm ~ 0.1 -0.5 mJy. At longer 
wavelengths, a large part of the CIB is due to 
faint 24 /im sources (624^™ °° 0.06-0.3 mJy) as 
shown b y cross-correlat i ons w ith BLAST data at 
350 /im (jMarsden et all 120091) and SCUBA data 
at 450 and 850 /im ([Serjeant et~ai1l2008l ). Going 
from far-IR to sub-millimetric, the CIB becomes 
progressively dominated by more luminous and 



higher redshift objects (jLagache, Puget fc Dole 



20051 ). Indeed, optical spectroscopy have shown 



that more than 70% of these sub-mJy 24 /im 
sources which are major contributors to the CIB 
lie beyond z of 0.6, with a primary peak at 
z ~ 1 and a secondary one at z ~ 2, consis- 
tently with the modeling of d e ep 24 /im num- 
ber co unts ( Caputi et al. l2007t Le Floc'h et al. 



20051 : lLagache et al.1 120041 : iDesai et al.l 12007). 



These sub-mJy 24 /im sources correspond to lumi- 
nous infrared galaxies (LIRGs, ^8-iooo/xm °° 10 11 - 
1O 12 L ) at z~l and to ultra-luminous infrared 
galaxies (ULIRGs, ^8-iooo/xm ~ 1O 12 -1O 13 L ) 
at z~2. The evolution of the IR popula- 
tion at z~l-2 is clearly illustrated by the in- 
frared luminosity functions at these early epochs 



irared luminosity junctions at tnese early epocns 
( Le Floc'h et a l. 2005; Caputi et alJ2QQ7tlReddy et al 



20081 : iMagnelli et al.l 120091 ) In essence, the knee 
of the luminosity function, which marks the dom- 
inant contributors to the energy density, moves 
to higher luminosity with redshift: from Ljr ~ 
(3-5)lO n L at z~l to (3-5)lO 12 L at z~ 
2. As shown by p ublished studies (see, e.g. 
iLe Floc'h et al.l (|2005| )). z~l LIRGs are responsi- 
ble for 70 ± 15% of the total (UV+IR) luminosity 
density, and, at z > 1, the contribution of ULIRGs 
becomes increasingly dominant. 

While deep 24 /im surveys have revealed LIRGs 
and ULIRGs as primary contributors to the global 



energy production at early epochs, the fraction 
of their bolometric luminosity contaminated by 
AGN emission is not well determined. Stan- 
dard multi-wavelength criteria can be used (us- 
ing X-ray emission, optical spectroscopy, IRAC 
photometry) but none of them is able to deter- 
mine the relative AGN and starburst contribu- 
tions to the bolometric luminosity of IR-luminous 
galaxies (that is mostly equal to the 3-1000/im 
luminosity). As of today, a very efficient way 
to probe dust obscured AGN emission, is via 
deep mid-IR spectroscopy. The difference be- 
tween the mid-IR spectra of a starburst and an 
AGN is striking. Starbursts are often character- 
ized by strong, low-excitation, fine-structure lines, 
prominent PAH features and a weak 10 /im con- 
tinuum whereas AGN display weak or no PAH 
features, p lus a stron g mid-infrared continuum 
(see, e.g., |wt| (119921 )). Mid-IR spectroscopic 
surveys have thus been conducted on a number 
of low to high redshift sampl es of IR-luminous 
galaxies with th e Spitzer/IRS (lArmus et al. 20071; 
Yan et al.1 [20071 : ISajina et all 12 0081: iFarrah et al ' 



20081: IValiante et all [20071: ISpoon et all 12007 



Wee dman et al. 112 0061 : Murphy et al .1120091 : 1 Veilleux et al 



20091 ). However, these samples are limited to the 
bright-end of the IR luminosity function. Previous 
deep mid-IR spectroscopic observations have been 
attempted at the faint 24 /im flux level of 0.1- 
0.5 mJy, bu t the number of observed sources is 



very small (|Teplitz et al.l 120071 : iBertincourt et al. 



l2QQ9h . To directly determine the AGN fraction at 
the lowest IR luminosity possible, we carried out 
an observational program to obtain low resolu- 
tion, mid-IR spectra of 48 galaxies — 24 at z ~ 1 
and 24 at z ~ 2 — with 24 /im flux densities ~ 0.1 - 
0.5 mJy. The main scientific goal of this study is to 
provide the first observational data to characterize 
the spectral properties of LIRGs and ULIRGs that 
are responsible for the bulk on the CIB and output 
energy at high redshift s {z ~ 0.7-2.5). Combining 
our deep mid-IR spectra with the available opti- 
cal and near-IR data in CDFS, we compare global 
properties such as colors, morphologies, and stellar 
masses of our sample galaxies with those of op- 
tically selected galaxies to provide observational 
constraints on the mid-IR properties of optically 
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selected IR faint galaxies at z~l-2. 

Throughout the paper, we adopt ^m = 0.27, 
= 0.73 and Ho = 71 kms -1 Mpc -1 cosmology. 

2. Data Sample 

2.1. Sample Selection 

LIRGs and ULIRGs are responsible for most 
of the integrated infrared background. To study 
these two populations of galaxies, we selected two 
samples of 24/im sources fainter than 0.5 mJy in 
two redshift bins (0.76 - 1.05 and 1.75 - 2.2). In 
these two redshift bins, the 24/im fluxes trans- 
late to infrared luminosities roughly in the ranges 
for LIRGs at z~l and ULIRGs at z-2. This 
sample is therefore crudely luminosity selected, 
and we do not apply any other selections. The 
flux lower limit (~ 0.2 mJy) is dictated by the 
sensitiv ity of the Spitzer InfraRed Spectrograph 
(IRS) ([Houck et all 120041 ). as demonstrated by 



the IRS instrument support team (jTeplitz et al.l 
120071 ). These targets were selected mostly using; 
photo metric redshifts computed by ICaputi et al ' 



(|2QQ6h since only 50% and 10% of them had 
spectroscopic- z in the z~l and 2 bins, respec- 
tively. We note that all sources but two are actu- 
ally in the chosen redshift bins, as confirmed by 
their mid-IR spectroscopic redshifts (see § 14. 2p . 
Our 48 galaxies were selected in the Chandra 
Deep Field South (CDFS) for which a large sam- 
ple of public ancillary data is available. The se- 
lected ULIRGs are all the 24 /im sources in the 
Chandra Deep Field South GOODS area with 
S(24fim)~ 0.14-0.55 mJy at z = 1.75-2.4, ex- 
cluding three sources too close to a neighbor 
source to obtain unblended spectra. The LIRGs 
have been selected among 24/im sources with 
S(24/ira)~0.2-0.5mJy at z = 0.76 -1.05 in order 
to avoid a large- scale structure localized at z 0.734 
( Szokoly et al.ll20o3 ). To have a more represen- 
tative sample of LIRGs, spanning different lumi- 
nosities within the LIRG range, we included in our 
selection sources from the "extended CDFS" area. 
The spatial distribution of the targets relative to 
the available ACS images is shown in Figure [U 
Tableland Table [3] list the 48 targets with J2000 
coordinates, on-target integration times for vari- 
ous slit modules, 8, 16, 24 and 70 /im flux densities 
and redshifts. 



One important question, which is also the 
basic driver of this program, is how represen- 
tative our sample is at the two redshift slices. 
Figure [2] quantitatively illustrates that our tar- 
gets represent populations of infrared luminous 
galaxies which are energetically important in 
comparison with other sources. We plot the 
contribution of 24 /im sources to the co-moving 
monochromatic 12 /im and 8 /im luminosity den- 
sity per logarithmic luminosity interval as a func- 
tion of luminosity at z~l and 2 respectively. 
Each of our targets is marked on the luminos- 
ity scale with a dashed line. The integrated co- 
moving lumin osity density is comp uted using both 
model (line) ([Lagache et al.ll2004l ) and previously 
publis hed, measured infrared luminosity func- 
tions (ICaputi et al.l 120071: iLe Floc'h et al.l 120051 : 
iMagnelli et al.l 120091 : iReddy et al.ll2008l ). It is in- 
teresting to note that at z ~ 1, the combined lu- 
minosity functions independently derived by three 
different groups suggest that the infrared luminos- 
i ty density turns over around L^lmp- ^ 10 10 " 7 ^ 



(ILe Floc'h et al .11200 d ICaputi et alJ2QQ7HMagnelli et al 



20091 ). Although the knowledge of the IR lumi- 
nosity function is still poor at z ~ 2, particularly 
at low luminosities, the critical luminosity range 
Lg/xm ~ 1O 11_11,5 L0 is fairly clearly highlighted in 
in the right panel of Figure [2] We conclude that 
our targets indeed probe the infrared populations 
which dominate the energy productions at these 
two cosmic epochs. 

Except having higher dust obscuration, our 
sample galaxies do not have extremely deviant 
properties in the rest-frame UV/optical in com- 
parison with galaxies selected at observed optical/near- 
IR band. For example, § 14.71 discussed in detail 
that the observed optical/near-IR colors of our 
sources are very similar to those of extremely red 
galaxy populations selected by large area K-band 
surveys. Figure [3] shows the most recently pub- 
lished stellar mass functions for redshift slices of 
0.8-1.0 (black lines and circles) and 1.5 -2.0 (red 



lines and circles) bv lllbert et al.l ((2009), and the 
dashed, vertical lines mark the estimated stellar 
masses for the galaxies in our sample. The method 
used to derive the stellar masses is described in 
It assume s a Kroupa's initial 



Wuvts et all (|2 




mass function (jKroupa 2001), stellar population 
synthesis models by iBruzual fc Chariot I (2003) 
and the reddening law of ICalzetti et al. ([2000). 
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Fig. 1. — IRS targets overlapped on ACS images from GOODS (light gray shaded) and GEMS (gray shaded). Empty and 
full circles mark the z ~ 1 and z ~ 2 samples, respectively. The two internal contours mark the border of t he 16^m imagin g 
(deep and ultra-deep observations). The external solid contour shows the coverage of deep Chandra imaging (|Luo e t al. 2008), 
while the external dashed contour marks the field covered by COMB017 ( Wolf et al. 200J) . 



The same set of assumptions are also used by 
Ilbert et al.l ([2009). This figure illustrates clearly 
that our z ~ 1 LIRGs mostly have stellar masses 
near or a bit lower than M*, the turn-over of 
the mass function, thus probing the most repre- 
sentative and important mass scales at this red- 
shift slice. The z ~ 2 ULIRGs reaches much lower 
masses, < 10 n M^ , than what was published in 
Ilbert et al.1 (|2QQ9h using COSMOS data. This is 
because the complete catalog used to construct 
these mass functions is cut at %ab ^ 25.5. This 
suggests that many of our z~2 ULIRGs have 
fainter optical magnitudes {%ab > 25.5). 



2.2. IRS Observations and Ancillary Data 

As shown in Figure HJ our targets have 24 /im 
fluxes of 0.14-0.55 mJy in comparison with the 
distribution of all 24 /im galaxies within the 
GOODS-South field (approximately 0.05 sq degs). 
For these ultra-deep, low resolution observations, 
we use the IRS in mapping mode, with 4 dithered 
positions along the slit for each target. All ob- 
servations were scheduled during each campaign 
just after the dark calibration observations. This 
method has been shown to be effective in reducing 
potential latencies left by previous o bservations 



of bright sources dTeplitz et al.l 120071 ). Our IRS 



spectra will sample important PAH features at 
6.2 /im and 7.7 /im, with the observed coverage of 
7-21 /im for LIRGs and the 14-35/im range for 
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Fig. 2. — Co-moving 12 and energy output per dex of 12 and 8/im luminosity (L dN/dLog L) at z = 0.94 and 2 = 1.91, 

corresponding to the mean reds hifts of our two bins. Continuous line s are fr o m the model of Lagach e et al.l (|2004h . Diamonds, 
squares, and triangles are from lLe Floc'h et al.l (120051 ): [Caput i et al ] j2007MReddv et al.l (|2008l), respectively. At zy 1, stars 
and crosses are derived using the 0.7 < z < 1 and 1 < z < 1.3 15 jum luminosity functions from Magnelli et al. ( 200§). Dotted 
lines represent the rest-frame 12 and 8 fim luminosities of our 48 targets computed from their 24 fim fluxes. No color correction 
has been applied. 




Log lo [M star (M )] Log 10 [M star (M Q )] 

Fig. 3. — Stellar mass functions at z~ 1-2 published by[nbert et al. (2009]using the COSMOS data. The dashed lines mark 
the estimated stellar masses of our sources estimated as in IWuvt s et al. (200|). 



Fig. 4. — 24 /im fluxes of our IRS targets (shaded his- 
togram) compared with the flux distribution of all the 
24 /im sources in the GOODS-CDFS field (solid histogram). 



ULIRGs. 

The CDFS field is remarkable for the great 
amount of ancillary data publicly available. In 
this paper, we make use of this rich dataset, in- 
cluding Chandra X-ray data, deep HST images 
at optical/UV, ground-based near-IR data, and 
the Spitzer images taken with MIPS (24 /im and 
70 /im), IRS (16 /im), and IRAC (8.0 /im). We 
reduced the Spitzer data by ourselves, and di- 
rectly to ok the reduced HST images from the 
GOO DS dGiavalisco et al.1 l2004h and GEMS sur- 
veys ( Rix et al.112004 ), as we ll as the recent FIRE- 
WORKS ( W uvts et al.ll2008h and MUSYCflmulti 



wave length photometric catalogs (Ga wiser et al 



20061). Most of the sources of our IRS sample are 
covered by the 2 Ms Ch andra Deep-Field South 
Survey (|Luo et aD l2QQ8h fsee Figured])- We use 
the X-ray catalog from iLuo et al.l (j2QQ8h for our 
analysis in 



2.3. Optical Spectra 

Published optical spectra exist for 17 ULIRGs 
and 9 LIRGs of our sample. The redshift estimate 
is usually easy for the z ~ 1 sample because of the 



large number of features in the spectra ([OIII], Ca, 
G-band) and becomes more complex in the case of 
z ~ 2 sources. In fact, there are good optical spec- 
z estimates for 9 LIRGs while only 7 ULIRGs have 
optical spectra with enough features to give reli- 
able redshifts (see Table [3]). Remarkably, all these 
spectra are available as FITS files and we were able 
to inspect each of them and confirm or correct the 
redshift estimates of the original publications. In 
the case of L4900, th e published z spec (opt) of 0.978 
(Popesso et al . 2008) was revised to 1.045 which is 
now perfectly consistent with the estimate derived 
from the IRS spectrum. All the spectra with good 
redshift estimates (according to the published ref- 
erences) are reported in Figure T he only excep- 



tion is the spectrum of U5795 (from lMignoli et aL 



(|2005h ) whose FITS file is corrupted. The pub- 
lished optical redshift is not reported in Table [3] 
since we were not able to confirm it. On the other 
hand, we include the spectrum of U4631 in the 
figure and table although it has been classified as 
unsure since it has clear Mgll and Fell lines and 
its optical spec- z is similar to that of estimated 
from the mid-IR spectrum. The optical spectral 
lines will be used in the discussion on the AGN 
fraction in § 14.31 

3. Mid-IR Data Reduction 

The reduction of the IRS data was carried out 
with a new software written by the first author 
specifically for low resolution IRS spectra. A de- 
tailed description of the technique used is given in 
the appendix. 

We reduced and measured fluxes of our targets 
in the 8 /im IRAC, 16 /im IRS, 24 and 70 /im MIPS 
images. For the IRAC data, we started from the 
Spitzer BCDs applying artifact corrections before 
mosaicking. The artifact corrections include col- 
umn pull-downs, muxbleed, optical banding, and 
background droops. We mosaicked the image with 
MOPEX with a pixel size of 0.6". The photome- 
try within an aperture o f 3.8 in radius was mea - 
sured using SExtractor ( Bert in fc Arnoutslll996f ). 
The final IRAC 8 jam flux is computing using a 
flux conversion factor of 0.2021 MJy/sr/(DN/s), 
and aperture corrections of 1.84 (see SWIRE tech- 
nical reporld, page 31). For the IRS 16 /im im- 
ages, we obtained a superflat by stacking the orig- 



1 www. strw.leidenuniv.nl/fireworks, www.yale.edu/MUSYC 



2 swire . ipac . caltech . edu / swire / astronomers / publications 
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Fig. 5. — Optical spectra available for LIRGs (top) and ULIRGs (bottom) of our sample. The main absorption and emission 
features are marked. Estimated optical redshifts are reported between brackets. 
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A/(l + z) [jum] 



Fig. 6. — Observed spectra for our sample of 24 LIRGs. The solid line, corresponding to the data with a smoothing of 5 
pixels, is on the top of the shaded la err or bars for each pixe l. Broad band fluxes are overplotted. Do tted lines are either 
the averaged local starburst template from Brand fet al.l (|2006h or the hot dust dominated spectrum from lSpoon et all (|2007l ) 
rescaled to the data. Source names and redshifts appear on the of each panel. 



inal BCDs. Each BCD is then corrected with tracted the s ources using PSF fitt ing with the code 

this superflat. Because we are interested in only StarFinder (|Diolaiti et al.l |2000[ ) inside an aper- 

point sources, we subtracted the median back- ture of 9.45 arcsec and applied an aperture correc- 

ground from each BCD before mosaicking them tion of 1.1626 (computed with a theoretical PSF 
with MOPEX with a 0.9 arcsec pixel. We ex- 
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Fig. 6. — continued. 



obtained with STinyTinU). The flux conversion 
used is 0.0117MJy/sr/(DN/s). In the case no 
counterpart was found, we computed also 5a flux 
upper limits. 

The 24 /im image was reduced usin g an im- 



proved version of the method described in Fadd a et al 



3 ssc.spitzer.caltech.edu/archanaly/contributed/stinytim 



(2006), which now corrects the long-term tran- 
sient and the regular short-term baseline varia- 
tions. The long-term is treated as a multiplica- 
tive effect since it has been proved in the case of 
MIPS 70 and 160 /im for which we have an inter- 
nal calibration that it is due to a variation in the 
responsivity of the detector. The short-term reg- 
ular variation of the baseline is due to distortions 
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Fig. 7. — Observed spectra for our sample of 24 ULIRGs. The solid line, corresponding to the data with a smoothing of 
5 pixels, is on the top of the shaded la er ror bars for each p ixel. Broad band fluxes are overplotted. Do tted lines are either 
the averaged local starburst template from lBrandl et a l. (2006) or the hot dust dominated spectrum from[Spoon et al.l ( 2007) 
rescaled to the data. Source names and redshifts appear on the top of each panel. 



introduced by the different positions of the cryo- 
genic mirror assumed to freeze the image during 
a scan. This variation is corrected by taking into 
account these distortions. We made the mosaic 
image with MOPEX and extracted sources using 



StarFinder with an aperture of 14 , and apply- 
ing the aperture an d color correction of 1.158 and 
0.961 respectively dFadda et al.ll200 6lFl. The flux 



4 The MIPS photometric system is tied to spectra with con- 
stant vfv, see MIPS data handbook, page 31 



10 



0.5 



U4812 [1.910] 





U4950 [2.291] 



.ImmImmII 



U5150 [1.898] 




0.5 - 



1-3 



0.5 



U5152 [1.794] 



U5795 [1.703] 



U5632 [1.998] 



U5652 [1.616] 






U5801 [1.843] 



U5805 [2.073] 






TT 



0.5 



U5829 [1.742] 



U5877 [1.886] 





U16526 [1.749] 




5 6 7 8 9 10 



6 7 8 9 10 

A/(l + z) [/am] 



6 7 



10 



Fig. 7. — continued. 



conversion is 0.0454 MJy/sr/(DN/s). 

Finally, the 70 fim imaging data is the deepest 
available, from the Far-Infrared Deep Extragalac- 
tic Legacy Survey (FIDEL jf). We start our reduc- 
tion from the raw data using the off-line Spitzer 



5 http:/ /ssc. spitzer. caltech.edu/legacy/all. html 



software (GeRI0). The pipeline corrects the cali- 
bration drifts, non-linearity, and latencies due to 
stimflashes, and applies the flux conversion before 
computing the slopes of the ramps. We remove 
any streaking in the final images by using filter- 
ing on the first-pass images with bright sources 



7 http: //ssc. spitzer. caltech.edu/mips/gert 
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Fig. 8. — Available broad-band fluxes and IRS spectra for three LIRGs without strong PAH featu res. The photo metric 
redshifts used in the sample selection were 0.754 and 0.776 for L5511 and L6211, respectively (|Caputi et a l. 2006). The 
proposed redshifts of 2.08 and 1.84 for L5511 and L6211, respectively, correctly place the valley between the parts of the SED 
dominated by stellar and dust emission as well as a strong feature visible in the L6211 IRS spectrum (identified as 6.2 /im PAH). 



masked out. The final mosaic with a pixel size of 
4 was obtained with MOPEX and sources were 
extracted with StarFinder with an aperture of 42 
in diameter. We apply the aperture correction of 
1.1665 (computed with STinyTim) and a color cor- 



rection of 0.918 to tie the photometric system to 
spectra with constant vf v . The flux conversion is 
702MJy/sr/(DN/s). In the case no counterpart 
was found, we computed also 5cr flux upper limits. 
The fluxes in these bands are reported in Table [3] 
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and in Figures [6j [7J and [T3J The flux limits in 
these two figures are in 5<r for all bands. All the 
flux conversion factors reported are the conversion 
factors applied by the Spitzer pipeline used at the 
time of our data reduction. We report these fac- 
tors, as well as the aperture corrections used, to 
allow a comparison with older or future reduction 
of other Spitzer data. 

4. Results 

4.1. Observed Mid-IR Spectra 

Figures [6] and [7] present the observed mid-IR 
spectra for the 24 LIRGs at z ~ 1 and 24 ULIRGs 
at z ~ 2. In each panel, the smoothed spectrum in 
the rest frame wavelength is shown as a red solid 
line, a local spectral template as a dashed line, and 
the unsmoothed spectrum with errorbars as an or- 
ange histogram. Each panel indicates the broad 
band fluxes at 8 /im, 16 /im, and 24 /im, the object 
name, the redshift based on the mid-IR spectrum, 
or optical redshifts when the mid-IR spectra do 
not yield redshifts (which is only the case for 3 
LIRGs). Downward arrows mark 5<r flux upper 
limits. The spectral template shown in dashed 
line in each panel is ei ther the aver aged starburst 
spectrum obtained by iBrandl et al.l (j2QQ6h based 
on the Spitzer spectra of 22 l ocal starburst gala x- 
ies, or the AGN template 1A (|Spoon et al.ll2QQ7h if 
the observed spectrum has no obvious PAH emis- 
sion. Redshifts are estimated from the mid-IR 
spectra by cross-correlating these two templates 
to the data. On Figures [6] and [7J the overlapped 
local spectral templates allow a simple characteri- 
zation of the observed spectra of our high-z LIRGs 
and ULIRGs. 

4.2. Mid-IR Spectroscopic Redshifts 

Of our 48 sources, we were able to estimate 
redshifts for 44 of them from their mid-IR spec- 
tra based on the PAH emission features and sil- 
icate absorption. We estimated the redshifts by 
cr oss-correlating the a veraged starburst template 
of IBrandl et~alT(|2QQ6h . For three sources in the 



1.2 



LIRG sample, we were not able to estimate spec- 
troscopic redshifts from the IRS data because of 
the lack of recognizable features. When plotting 
their spectra with the broad band photometry 
points (see Figured]), it is clear that two of them 
(L5511 and L6211) are at redshifts higher than 
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Fig. 9. — Comparison between mid-IR spectroscopic 
redshifts and optical redshifts (spectroscopic and photo- 
metric). In the top panel, photometr ic redshifts are from 
IWuvts et all ([20081 ) (full circles) and ICaputi et all (|2006h 
(empty circles). 



those estimated bv lCaputi etHI (|2QQ6h (0.754 and 
0.776, respectively). They are actually ULIRGs at 
z = 2.08 and 1.84, respectively. Furthermore, Fig- 
ure [6] shows that the spectrum of L6211 has fairly 
strong 6.2 /im PAH emission at z = 1.84, with the 
rest of the PAH spectrum not sampled by the ob- 
servations (LIRG observations used only LL1, not 
LL2). The wrong estimates are due to wrong iden- 
tifications with optical sources at ~ 1.5 arcsec from 
their positions. The real counterparts are faint 
sources very close to these bright galaxies and usu- 
ally blended with them in ground images which 
are not reported in the COMBO 17 catalog. The 
faint optical counterparts appear very clearly in 
the HST images (see Figure [24]). 

Figure [9] compares mid-IR spectroscopic red- 
shifts with two sets of photometric redshifts in the 
top panel and with the optical spectroscopic red- 
shifts in the bottom panel. The mid-IR spectro- 
scopic redshifts are consistent with the available 
optical spectroscopic and photometric redshifts 
within 20%, and the photometric redshifts have 
larger uncertainties than that of mid-IR spectro- 
scopic redshi f ts. T he photometric redshifts from 
Wuvts et al. (|2008h have a slightly lower disper- 
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sion than those from ICaputi et al.l (|2QQ4f ) since 
they are based on optical and infrared data. Ta- 
ble [3] reports preferentially photometric redshifts 
from the former study. We conclude that our sam- 
ple selection based on initial photometric redshifts 
has proven effective. 

Our IRS spectra have provided 8 n e w spe ctro- 
scopic redshifts used in Wu vts et al.l (2009) for 
which ground-based optical/near-IR spectroscopy 
is difficult. Particularly, at z~l. 7-2.0, H a falls 
between the H- and K-band atmospheric windows, 
making ground-based near-IR spectroscopy im- 
possible. For 24 /am galaxies in this redshift range, 
their Spitzer mid-IR spectra provide useful dis- 
tance measurements. 
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4.3. AGN fraction 
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Fig. 10. — Hard X-ray 2-8 keV luminosity versus equiv- 
alent width of the 6.2 fim. PAH feature for the IRS sources 
with X-ray counterparts. The infrared output of the ob- 
jects on t he left of the vertical line is dominated by AGN 
emission (|Armus et al . 2007; Veilleux et al. 2009). Objects 
above the horizo ntal line are t ypica lly AGN dominated 
sources (see, e.g., Fadda et al. (2002)). The dashed ver- 
tical line marks EWQ.2^m = 0.2 |Um, below which galaxies 
are AGN dominated. See § 14.31 for the detailed discussion 
of this value. 



One motivation for this program is to use mid- 
IR spectral diagnostics to determine the contri- 
bution of AGN to the total energy output among 
LIRGs at z ~ 1 and ULIRGs at z ~ 2, particularly 



Fig. 11. — 6.2 /im equivalent width versus redshift. The 
lower limits come from upper limits to the continuum from 
16 /im observations. The horizontal line corresponds to 
the limit between AGN and star formation dominated sys- 
te ms. EWQ.2um = 0.2 fim i s derived from the combination 
of lArmus etahl (120071 ) and lVeilleux et al.l (120091 ). 



the obscured AGN which may not be detected by 
X-rays. It is important to understand that, al- 
though in some cases AGN can completely dom- 
inate the bolometric luminosity, there is an en- 
tire population of galaxies with enshrouded AGN 
whose bolometric luminosities are primarily from 
star formation. To objectively classify AGN, we 
use several indicators such as: broad lines and high 
ionization lines in optical spectra, lack of a 1.6 jam 
stellar bump in the SED, X-ray bright sources, 
low mid-IR 6.2 /am PAH equivalent width, and fi- 
nally their optical morphology. We note that our 
samples at z ~ 1 and z ~ 2 now have 22 and 26 
sources respectively, because two galaxies initially 
assumed to be LIRGs at z ~ 1 have now been de- 
termined to be ULIRGs at z~2 (see §E2J). 

The 2 Ms Chandra X-ray survey covers all the 
ULIRGs and most of the LIRGs of our sample. 
Six of our sour ces have X - ray co unterparts from 



the catalog of iLuo et al.l (|2QQ8f ) (see Table [Q) 



Two of them, U4950 and U5877, are clearly de- 
tected in the (2 - 8) keV hard X-ray band. An- 
other one (U4642) is detected at the limit of sen- 
sitivity of the survey. Two more ULIRGs (U5775 
and U428) are detected only in the soft X-ray and 
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one LIRG (L5420) is detected at a low-a level 
in the soft X-ray (secondary catalog). A cross- 
m atch with the more e xtended X-ray observations 
of iLehmer et al.l ([20051 ) does not yield any further 
counterpart. Figure [12] compiles all of the broad 
band photometry and mid-IR spectra for these 6 
sources. U4950 and U5877 have IRS spectra com- 
patible with a power law and also their broad band 
photometry is typical of an AGN, i.e. no 1.6 /im 
stellar bump and a power-law spectrum over the 
entire wavelength range. We note that the HST 
morphologies of U4950 and U5877 show compact, 
unresolved nuclei, supporting the AGN classifica- 
tion from both X-ray and mid-IR spectra. The 
other source detected in the hard X-ray (U4642) 
has, on the contrary, a more complicated morphol- 
ogy with multiple clumps. So, one clump can host 
the AGN, but the AGN does not dominate the 
optical-IR emission. All the other sources (U5775, 
U428 and L5420) have a softer X-ray emission, the 
mid-IR spectra have obvious PAH emission (see 
the quantitative numbers below) and their broad- 
band SEDs shown in Figure [12] reveal a strong 
stellar bump at rest-frame wavelength 1.6 /im, sug- 
gesting that although these three sources may have 
low-luminosity AGN, their host galaxies dominate 
the stellar emission in optical to near-IR. The X- 
ray luminosity of the fainter X-ray sources is also 
compatibl e with star formation using t he rel ation- 
ships from lRanalli. Comastri fc" Settil ([20031 ). 

The optical spectrum of U4958 (see Figure [5]) 
shows a broad CIII line and strong NV and CIV 
emission lines. Considering that the infrared spec- 
trum is rather featureless, we can classify this 
source as AGN dominated. 

The strength of PAH features has been widely 
used as a star format i on indicator (|Lutz et al 



used the code PAHFIT ([Smith et al.l 120071 ) to fit 
lines and dust features together with a continuum 
from dust and stellar emission to our IRS spectra. 
During our analysis, we discovered two bugs in 
the routine that computes equivalent widths. The 
routine was incorrectly applying Gaussian profiles 
instead of Drude profiles in the computation and, 
more importantly, using the flux of the line cor- 
rected for extinction and not corrected continuum. 
These bugs have been reported to the author of 
the code and will be fixed for the next release. To 
measure the equivalent width of the 6.2 /im PAH, 
we used a Chiar-Tielens silicate profile and added 
an absorption component to take into account 
the water ice and hydrocarbon f eature s at 5.7- 



7.8 /im. Following IVeilleux et al.l ((2009), we used 



19961: IGenzel et afll998t lArmus et alJ2QQ7HVeilleux et al, 
2009h . Particularly, the 6.2 /im PAH is relatively 
isolated, and its rest-frame equivalent width can 
be measured cleanly and used to classify starburst- 
dominated systems. This is supported by the fact 
that EWe,2fim is shown to broadly anti-correlate 
with flux ratios of [NeVl/fNelll and [OVl/fNelll 
( Lutz et_al[ll996l: IGenzel et al.lll998t lArmus et at 



20061 : IVeilleux et al.l 120091 ). where the detections 
of [Ne V]A14.3/im and [O IV]A25.9/im indicate 
the excitation of the ionized gas by black hole 
accretion due to their extremely high ionization 
potentials (97.1 eV and 59.4 eV respectively). We 



the profile taken fro m observations of F00183-7111 
([Spoon et al.l |200J) for this absorption feature. 
The addition of the water ice and hydrocarbon 
absorption features at 5.7- 7.8 /im is important 
mainly in strongly absorbed local ULIRGs (see 
§ 14. 6|) . Neglecting this component in our sample 
spectra and other high-z spectra does not change 
significantly the EWQ.2/j,m estimates. 

In the case of LIRGs, the signal-to-noise ratio of 
the spectra is high enough to measure directly the 
equivalent width of the 6.2 /im PAH feature. In 
the case of ULIRGs, the low continuum combined 
with the higher noise of the spectra lead to large 
errors in the estimates of EWe^fim- In particu- 
lar, the error is dominated by the uncertainty in 
the estimate of the continuum from the spectrum. 
To obtain firmer estimates for individual sources, 
we used 16 /im fluxes as estimates of the contin- 
uum since they have a much lower uncertainty. 
We have 16 jam fluxes for 14 ULIRGs and upper 
limits for an additional 9. Although, in general, 
most of the 16 jam flux is due to the continuum, 
the contribution of the 6.2 /am PAH is increasingly 
important at redshifts lower than 1.7. To com- 
pensate for this effect, we applied a multiplicative 
correction to the 16 jam flux equal to the ratio of 
the 5.8 jam flux and the expected 16 jam of the av- 
eraged starburst spectrum moved at the redshift 
of the ULIRG. This correction ranges from 0.6 at 
to 0.9 at z~2 and does not vary signifi- 
cantly by using different templates as the average 
HH-like ULIRG spectrum which better describes 
ULIRGs around the Si absorption. We were able 
to measure the EW^^^m for all the sources with 
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Table 1 

Cross-correlation with X-ray survey 



Name Distance Flux [10 16 erg /s/ cm 2 ] L2-8keV 

(arcsec) (0.5-8 keV) (0.5-2 keV) (2-8 keV) [10 43 erg/s] 



primary 


catalog: 










U428 


0.29 


0.76 


0.35 


<1.26 


<0.27 


U4642 


0.12 


2.05 


0.74 


1.41 


0.35 


U4950 


0.18 


60.90 


13.00 


47.00 


18.9 


U5775 


0.48 


<0.97 


0.42 


<2.39 


<0.35 


U5877 


0.18 


36.40 


0.73 


36.40 


9.1 


secondary catalog: 










L5420 


0.40 


<2.07 


0.63 


<2.91 


<0.18 



the exception of U4958 which we already classified 
as AGN on the basis of its optical spectrum. 

The 6.2 /am equivalent widths of the X-ray emit- 
ting sources are reported in Figure [T0l and support 
our previous discussion. Studies of local LIRG 
and ULIRG spectra have shown that the per- 
centage of AGN contribution to the total energy 
outpu t correlates with EW^^^m- lArmus et al. 
(|2007l ) estimates that E We> 2um < .3 /am are 
AGN-dominated, where Veilleux et al. (|2QQ9h put 
the limit around 0.1 /am while at 0.3 /am the 
AGN contribution is around 30%. The two re- 



sults are not in contradiction since Veilleux et al. 



(2009) compare AGN fractional contributions to 
the bolometric luminosities by applying bolomet- 
ric corrections to the ir numbers (see Table 10 in 



Veilleux et al. d2009l) and disc ussion in the Ap- 



pendix). lArmus et al.l ([20071 ) compare instead 
AGN fractional contributions to the [Nell] and 
mid-IR luminosities. Nevertheless, a ccord ing to 

also 



other criteria used by I Veilleux et al.l (|2C 
sources with EW^^nm up to 0.2 jam can be AGN- 
dominated (see Fig. [18]). In the following, we 
assume 0.2 /am as a conservative limit. 

Figure [10] shows that the mid-IR 6.2 jam PAH 
equivalent criteria has yielded the same conclusion 
as we stated above, that 2 X-ray sources U4950 
and U5877 are strong AGNs, and the rest of four 
have their infrared emission dominated by star for- 
mation. 

Since many LIRGs are outside of the X-ray field 
or in the less sensitive parts and because of the lim- 



ited hard X-ray sensitivity, it is useful to analyze 
the 6.2 /am equivalent widths of our entire sam- 
ple. Figure [TT] shows the 6.2 /am equivalent width 
versus redshift for all 48 sources with measurable 
6.2 /am PAH. The infrared emission of almost all 
the LIRGs and the majority of the ULIRGs is 
dominated by star formation. The LIRG with low 
6.2 /am equivalent width (L3945) lie in the less sen- 
sitive and external regions of the X-ray survey and 
therefore it is not a surprise that it is not detected. 

Among LIRGs, therefore only one source is 
dominated by the AGN emission (5% of the sam- 
ple). In the case of ULIRGs, we exclude L5511 
(reclassified as ULIRG) from the sample because 
its spectrum does not have enough coverage to 
yield usable 6.2 /am equivalent width measure- 
ment. Thus, of the total 25 ULIRGs at z ~ 2, only 
12% (3/25) have a dominant AGN. 

These numbers are significantly lower than 
those for bright (S u > 1 mJy) 24 /am samples in the 
same redshift range. Based on the mid-IR spectral 
features as well as IR SEDs, the fraction of AGN 
domin ating the IR lum i nositi e s are close to (50 - 
75)% (|Yan et al.l I2QQ5L 120071 : LSaiina et al.l I2007L 
Including weak AGNs or starburst/AGN 



composite systems, the majori ty 75%) of brigh t 
24 /am samples contain AGNs (|Sajina et al.ll2QQ8h . 
Studies of local ULIRGs have found that the AGN 
fraction becomes significant at Lir > 10 123 L© 
( Lutz et al. 1998), although the exact numbers 
suffer from many uncertainties. There are specula- 
tions (see § 14.6ft that at high-z, this AGN-starburst 
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Fig. 12. — Broad band photometry and IRS spectra for the six IRS sources with X-ray counterparts. Only two of them 
(U4950 and U5877) have SED and spectral features compatible with AGN dominated sources. 



transitional luminosity moves to higher values. 
One supporting evidence has been SMGs, whose 
Lir is on average > 1O 12,3 L , and it is commonly 
accepted that bolometric luminos ities of SMGs are 
mostly from starbursts (see, e.g., Alexander et al.l 
(j2QQ3h : I Valiant e et al.l ([20071 )). One result from 



this paper is that the AGN fraction for z ~ 2 
galaxies is small for Lir ~ 1O 12,6 L , which implies 
that at high-z, the AGN-starburst transitional lu- 
minosity must be higher than Ljr^ 1O 12,6 L . 
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Fig. 13. — Templates of average starburst (Bra ndl et all l2006) and ULIRG composites ( Vei lleux et afll2009D rescaled to the 
smoothed IRS spectra (red line) of non AGN-dominated ULIRGs. The best fit compatible with the 70 fim flux is shown in blue. 



4.4. Stacked Mid-IR Spectra 

How do the relative PAH strength of high-z 
LIRGs and ULIRGs compare with those of local 
infrared galaxies? As shown in Figures [6] and [7J 
our individual spectra are very similar to the av- 



eraged z ~ starburst spectrum ( Brandl et al.l 
2006). We note that our spectra are typically 
not deep enough to individually probe the silicate 
absorption feature at 9.7 /im. In this section, we 
make this comparison by producing stacked spec- 
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tra at z of 1 and 2 redshift bins. This technique 
has been widely used in recent years to detect and 
quantify certain weak features which are not de- 
tected in the data of individual sources. 

We exclude from the stacking the few sources 
whose emission is AGN-dominated: three ULIRGs 
(U4958, U4950, U5877) and one LIRG (L3945), 
as explained in § 14.31 We choose to normalize 
the spectra to their 24 jam rest-frame luminosity 
which correlates w i th th eir total infrared lumi- 
nosity ( Rieke et al. 2009f ). Here, the rest-frame 



24 /am luminosity is calculated differently for the 
two samples. For z ~ 1 LIRGs, we simply interpo- 
lated the 8, 24, and 70 /am fluxes in the log (A) ver- 
sus log(V F v ) plot. We do not consider the 16 /am 
flux since it is severely affected by the 7.7 /am PAH 



complex. In the case of z ~ 2 ULIRGs, the 70 /am 
flux is very close to the rest-frame 24 /am flux. Un- 
fortunately, for many ULIRGs we have only upper 
limits of the 70 /am fluxes. In these cases, we es- 
timate the rest-frame 24 /am fluxes by scaling a 
series of te mplates, including the averaged star- 
burst from iBrandl et al. (l2QQ6h and th e compos- 
ite ULIRGs from lVeilleux et al.l (|20Q9h . to fit the 
spectrum. The best fitting template compatible 
with the 70 jam flux or upper limit is then consid- 
ered to estimate the rest-frame 24 /am flux. The 
rescaled templates are shown in Figure [13l In Ta- 
ble OS we list the estimated total infrared luminosi- 
ties. 

The result of the stacking is shown in Figure [TU 
We have resampled the original spectra on a com- 
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Fig. 14. — The averaged z ~ 1 LIRG (left) and z ~ 2 ULIRG spectrum (right) obtained by stacking all the spectra which are 
not AGN-dominated. The spectra have been normalized to the 24 /im rest-frame luminosity before stacking them. The lower 
panel shows the number of spectra used for each wavelength. The continue and dashed green lines correspond, respectively, to 
the t otal flux and contin uum fitted with PAHFIT. The red and blue lines correspond , respectively, to the averaged starburst 
from iBrandl et all (120061 ) and the average local HII ULIRG from FVeilleux et al.l (|2009T ) normalized to the stacked spectra. 



mon wavelength grid as close as possible to the 
original data. At each wavelength, we computed 
the average values weighted by the errors on the 
single spectra. To check if outliers create artificial 
strong features, we also obtained spectra by using 
the biweight mean and we did find only negligible 
variations between the two methods. 

The averaged spectra show clearly the char- 
acteristics of starburst dominated galaxies with 
weak dust continuum at 4- 14 /am and strong PAH 
emission at 6.2,7.7,8.6, and 11.3 /am. Weak PAH 
features at 5.25 and 5.70 /am are also visible, es- 
pecially in the LIRG spectrum. For compari- 
son, we report in Figure [14] also local average 
spectra normalized to our data (by minimizing 
the x 2 °f the residuals). The average LIRG is 
surprisingl y close to the averag e starburst spec- 
trum from IBrandl et al.l (|2QQ6h . In the case of 
the average ULIRG, we co mpare it to th e aver - 
age starburst spectrum from IBrandl et all ( 20Q6h 



and to severa 


. averaf 


Veilleux et al. 


(120091 


Veilleux et al. 


(2009h 



The best fit among the 
spectra is obtained using the 
H Il-like ULIRG average spectrum. Although the 
average starburst spectrum fits well to the relative 



PAH strength in the averaged spectrum of z ~ 2 
ULIRG, the local average ULIRG reproduces bet- 
ter the 9.7 /am silicate absorption. This implies 
that high-z ULIRGs have more dust extinction 
than lower luminosity LIRGs. 

4.5. Comparison with Empirical IR SED 
Templates 

Empirical IR SED templates used to model 
counts from deep surveys and the cosmic infrared 
background are mostly based on local galaxies. 
One important question is how these empirical 
SEDs compare with our observed SEDs, includ- 
ing broad band photometry and spectra. Here we 
considered all the spectra which are not dominated 
by AGN emission and obtained composite spec- 
tra by weighting the individual spectra by their 
total infrared luminosities. The same procedure 
has been applied to the IRAC and MIPS fluxes of 
the sources. The composite SED plus IRS spec- 
tra are compared to four different models, among 
the most popular, in Figures [T5l and [T6l For each 
model, we chose the one with closest total infrared 
luminosity and rescaled it to the total infrared lu- 
minosity of the composite spectrum. In the case of 



20 



N 



10 



25 



10 



24 



10 



23 



10 



22 



10 



25 



10 



24 



10 



23 



10 



22 



- Dale & Helou 




: Chary & Elbaz 




_i i i i 1 1 1 1 



j i i i 1 1 1 1 



Lagache 




_i i i i 1 1 1 1 



j i i i 1 1 1 1 



_l I I I I I I I I I L_ 



Franceschini 




_l I I I I I I I I I L_ 



10 



100 1 
A[/xm] 



10 



100 



Fig. 15. — Compos i te LIRG spectrum and IRAC and MIPS fluxes compared to four popular theoretical templates 
(jFranceschini et alJlSool lLagache et al. 2004: lDale fcH elou 20021: ICharv fe Elbad l200lh . 



the Dale fe Heloul (l2QQ2h models, which give spec- 
tra as a function of the Fqq/Fiqq color, we used the 
relationship betwee n this color and the total in- 



frared luminosity in lChapin. Hughes fc Aretxaga 
(|2009h . 

Two models extend down to near-IR fluxes 

([Dale fc Heloul l2QQ2t ICharv fe Elbazll2QQlh. while 

the o ther two ([Lagache et al. 2004; Franceschini et al 
20091 ) are limited to wavelengths longer than 4 fim. 
It is interesting to note that of the four empiri- 
cal SED template models, only the Franceschini 
model templates fit the mid-IR spectral features, 



although all four models more or less in agree- 
ment with the overall observed SEDs for both 
LIRGs and ULIRGs. However, the few observa- 
tional constraints in the far-IR part are consistent 
with all four tem p lates. For LIRGs, only the 
Franceschini et al ] (|2QQ9h model SED template 



agrees with the observed spectral features and is 
able to reproduce the far-IR limit for the observed 
LIRG composite. 

Using the composite spectra plus the IRAC and 
MIPS fluxes, we tried also to estimate how much 
the presence of PAH features boosts the detection 
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Fig. 16. — Composite ULIRG spectrum and IRAC and MIPS fluxes compared to four popular theoretical templates 
dFranceschini et alJl2009l : lLagache et al. 2004: lDale fcH elou 20021: ICharv fe Elbadl200lh . 



of infrared galaxies in the 24 /am filter as a func- 
tion of redshift. As shown in Fig. [T7J we fitted the 
composite LIRG and ULIRG with PAHFIT to es- 
timate the continuum. Then, we computed the 
ratio of the flux in the 24 /im filter for the total 
spectrum and for the continuum only. The ratio 
of these two fluxes has two clear peaks approx- 
imately at z = 0.9 and 2 = 1.8 which corresponds 
to the redshifts where the 24 /im filter includes the 
two main PAH complexes (11.3 and 7.7/im). For 
comparison, the median redshifts of our two sam- 
ples are 0.95 and 1.87. Figure fT71 suggests that at 



z ~ 2, 24 /im broad band filter selection will prefer- 
entially select galaxies with strong PAH emission 
rather than power-law sources, and this bias can 
be as high as a factor of 2. 

4.6. 6.2/im equivalent width versus far-IR 
luminosity 

In this section, we explore the relationship be- 
tween two quantities which are mainly linked to 
star formation: the rest-frame 24 /im monochro- 
matic luminosity {vL v )2Ay,m an d the equivalent 
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Fig. 17. — Fraction of the total over continuum 24 /im flux for luminous galaxies in the range between redshift and 4, based 
on our composite LIRG and ULIRG spectra. The presence of PAH features boosts the detection of 24 fim sources around z ~ 0.9 
and z~ 1.8, the median redshifts of our two subsamples. The insets show the decomposition of the two composite spectra in 
PAH plus continuum and the locations of the 24 jum filter at the two sensitivity peaks. 



width of the 6.2 /im PAH feature. The rest-frame 
(vL v )2Aiim is related t o the total infrared luminos- 
ity Lir = L8-iooo//m ( Rieke et al.ll2QQ9h which is 
an indicator of star formation. On the other hand, 
PAH features are created in star-forming regions 
and destroyed around AG N, due to their intense 



i oniza tion fields ( see, e.g. JVoitl (|1992| ): iLutz et al. 



(1998)). We expect, therefore, that in presence of 
an AGN component, {vL^i^m is proportional to 
the sum of the AGN and starburst luminosities. 



On the other side, EW^^nm is the ratio between 
the strength of the line (due only to star-formation 
activity) and the underlying continuum (created 
by AGN and star- formation activity). These two 
quantities should behave like: 



24/im OC (LaGN + 

Lsb 
Lagn + Lsb 



(vL v ) 



(1) 

(2) 
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where Lagn and Lsb are the AGN and starburst 
luminosities, respectively. 

In Figure HH1 we report the 6.2 fim equivalent 
widths measured for the z~l LIRGs (blue dots) 
and ULIRG at z ~ 2 (black dots) as a function of 
their rest-frame 24 /im luminosities. We also plot 
the EWe^fim of the composite ULIRG (black dot 
with errorbar). The stacked ULIRG spectrum has 
a high enough signal-to-noise ratio to directly mea- 
sure its mid-IR continuum and therefore obtain an 
average EWe^m- 

In Figure fT8l we compare local starbursts 



sources have PAH emission stronger than that 
of local ULIRGs. This can be interpreted in at 
least two different ways. Either high-redshift 
sources are forming stars more efficiently than 
local counterparts or local ULIRGs have a higher 
AGN contribution to their infrared continua. Sup- 
porting evidence for the first hypothesis is that 
SMGs at z~2, which are the majority of the 
high-z luminous galaxies in our plot, have more 
gas, higher star formation rates and more efficient 



2001 



star formation than loca l ULIRGs (jGreve et al 



(Bran dl et al.M2Q06h , Uy ULIRGs (IVeilleux et al 



2009) with a collection of high-,? Spitzer ULIRGs 
and SMGs (Saiina et al. Menendez-Delmestre et al 



120071 : IPope et al.ll2QQ8t IValiante etal. 



JSooi 



We 



report also the average points from iDesai et al. 
(|2007l ) computed from local galaxies. To make this 
comparison meaningful, we re-measured EWe^^m 
and I/24 Mm for all of the sources included in Fig- 
ure [TBI We measured the EWe^m using our 
revised version of PAHFIT (see § !4.3p . The values 
of EW 6 . for highly extincted ULIRGs agree 



very w ell with those computed by IVeilleux et al 



(2009). On the contrary, for less extincted galax- 
ies we have systematically l ower values (~ 70%) 
probably due to the fact that lVeilleux et all ( 20091 ) 
do not fit directly the lines but uses templates to 
fit the spectra whose actual relative PAH strength 
may differ from that of templates. 

One striking result shown in Figure [18] is the 
maximum EW^^^m value, ~ 1 /im, an upper limit 
for all starburts, LIRGs and ULIRGs, indepen- 
dent of redshifts. This maximum limit is ex- 
pected from Equation [TJ indicating that when 
Lagn «Lsb, EW^^m should be constant. We 
have no explanation why this constant is ~1 /am. 
At L24 M m^lO n L , EW 6 . 2 ^m and L 2 4/xm are 
clearly anti-correlated with a slightly different 
slope for local and high-redshift sources. The 
relation is broad not only because of errors in 
the measurement of the EWe.2/j,m but also due to 
AGN contributions that dilute the EWe.2/j,m in- 
trinsically broadening the correlation. As a matter 
of fact, the dispersion of pure local starbursts (for 
which the AGN contamination is negligible) is 
much smaller than that of the local ULIRGs (that 
are contaminated by AGN). 

Among high luminosities sources {L2^ra > 
10 n L©), at a given luminosity, high-redshift 



20051 : iTacconi et alll2006l ). On the other hand, if 
the AGN contribution to the mid-IR continuum 
is smaller for z > 1 ULIRGs than that for local 
ULIRGs at a fixed luminosity, high-z ULIRGs 
would have bigger EWe^m due to less AGN 
boosting of the mid-IR continua. The good fit 
of the composed spectra of our sample with the 
averaged starburst template and the low estimated 
AGN contribution to their total emission (see also 
next section) tend to support this last hypothesis. 

4.7. Red Optical Colors 

Many of our sources have very red optical 
SEDs. As shown in Figure HH with definitions 
for Extremely Red Galaxies (ERG) and Dis- 
tant Red Galaxies (DRG) of (i - ij O Vegai >4.0 
and (J — iOveffa<^2.3, respectively ( JYan et al 



200< 



2005; iFranx et al.l 120031 : Ivan Dokkum et al 



2004D, we find that of the 24 z~ 2 ULIRGs, 63% 
(15/24) are ERGs and 33% (8/24) are DRGs. 
These near-IR color selected DRGs have very faint 
optical magnitudes, and their redshifts are diffi- 
cult to obtain wit h optical/near-IR spectroscopy 



(|Kriek et al.l 12008). Our IRS spectra have pro- 
vided 8 spectroscopic redshifts which do not 
have optical spec-z and have been used for the 
analys is of the re liability of DRG color selec- 
tions (|Wuvts et al.ll2009h . Our finding of 8 z~2 
ULIRGs having very red (J — K) colors suggest 
that a significant fraction of DRGs selected by Ri- 
band surveys are dusty star forming galaxies, not 
passively evolving stellar populations at z~2-3. 
This is consistent with previously published stud- 
ies, which found as much as 50 % of color selected 
DRGs being dusty starbursts ( Wuvts et al. 20071; 



2008; IWuvts et al. 



Labbe et alJl200a iKriek et al 
2009h . 

Furthermore, all of the z ~ 2 ULIRGs also sat- 
isfy the (z — B) vs. (K — z) color-color selec- 
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Fig. 18. — Equivalent width of the 6.2 /zm PAH feature versus the rest frame luminosity at 24 /imfor galaxies of our sample 
(2 ~ 1 and z rsj 2 galaxies are marked with blue and black dots, respectively) compared to local and other high-z galaxies. We 
note that we re-measure EWq^^tu and L<2^m for all objects in this figure. The bla ck dot with error bars is for the average z ~ 2 
ULIRG. The red triangles correspond to local starbursts from lBrandl et all (|2006f ), while the red circles are ULIR Gs from the 
local 1 Jy sample by Veill eux et al.l (|2009h . The empty circles are classified as at least 50% AGN dominated bv [Veilleux et al. 
(120091). Finally, t he green dots are high-red shift infrared sources from literature (Saiina et al 1 l20Q7l : Menendez-Delmestre et al. 
l2007l : lPope et al1l2008l : IValiante et al.ll2009l). All the equivale nt widths have been remeasured in the same way as the one from 
our sample. Cyan points and line are from lDesai et al.[(|2007T ). 



tions (BzK) (|Daddi et al.ll2004h . as shown in Fig- 
ure [20j These BzK color cuts shown as lines in Fig- 
ure [20] have been used to select star forming and 



passi vely evolved galaxies at z > 1.4 (jDaddi et al 
20041 ). These z~2 ULIRGs fall into the star form- 



ing BzK region on the (z-B) vs. (K-z) color-color 
diagram. As expected, their red BzK colors are 



mostly due to the dust absorption the UV pho- 
tons from O, B types of young stars produced by 
recent starbursts. 

In Figure [20j we report also the IR excess as 
a func tion of Ls^m computed as in iDaddi et al. 
(|2004h . All our z - 2 ULIRGs and a few LIRGs 



are IR excess galaxies according to the definition of 
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Fig. 19. — Color-color plots for galaxies in the GOODS field (dots) and galaxies of our sample. Full and empty circles mark 
the LIRGs and ULIRGs of our sample, respectively. 



Daddi et al. I (|2004h . Although it is a known fact 
that optical and IR surveys select galaxies with 
different amounts of dust, the question here is if 
our IR selected galaxies follow the same dust ex- 
tinction scaling relationship as optically selected 
objects with low dust content. Figure [21] shows 
the far-infrared to 1600A monochromatic lumi- 
nosity ratio as a function of UV spectral index 
/3, where f\~\P, 1000 A < A < 2500A. The ratio 
Lfir/L ° measures how much energy is re- 

160CL4 toJ 

distributed from UV to infrared, and j3 measures 
the UV SED slope, L.e UV reddening. Studies 
of local UV selected starbursts have shown that 
dust absorption is correla ted with UV reddening 
([Meurer et al.lll995l Il999h . This correlation pro- 
vides a powerful empirical tool, particularly for 
high-z UV/optical surveys, to recover the total, 
dust-absorption corrected, UV fluxes, using UV 
photometry alone. Figure [21] compares the lo- 
cati on of our sample with those of local galax- 
ies ([Meurer et al.l ll999) and high-redshift galaxies 



([Chapman et al 
plot L FIR /L 



20051 lR,eddv et al.l l2QQ6h in the 
versus j3. In the case of our 



160CL4 

sample, we used the available optical magnitudes 
from U to R to compute ^ 1600 ^ by interpolating in 
vF v , while was com puted using the relationships 
from IPapovich et al.1 (|2006h . The J^fir is com- 
puted from the rest-fr ame 24/im l u minos ity using 
the relationship from iRieke et al.l (l2009[) to com- 
pute the total IR luminosity and from lElbaz et al 
([20021 ). to pass from total IR luminosity to far- 
IR luminosity. In the case of SMG galaxies from 



Chapman et al.l ((2005), we used the available B 



and R magnitudes to compute 6 and L ° . The 

to 1 r 1600A 



I 1 J llUUU/l 

relationship from lPapovich et al.l ((2006) were used 
after applying median corrections to the B and R 
magnitudes to transform them in the appropriate 
bands. For the z~l sample, the UV luminos- 
ity and f3 are poorly constrained since we do not 
have UV observations. These sources appear to 
scatter between the normal star-forming galaxies 
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Fig. 20. — Left: optical/near-IR color plot; Right: IR excess vs rest-frame 8/im luminosity. In the two figures, sources from 
the ULIRG and LIRG samples are marked with filled and empty circles, respectively. The names indicate the two sources from 
the LIRG sample which are actually at redshift ~ 2. 



and the UV-selected starbursts relations, which 
as discussed above can be attributed to a range 
of UV extinction slopes. For our z ~ 2 ULIRGs, 
we note that they largely lie above the Meurer's 
relationship, although usually not quite so far as 
the SMGs and HLIRGs (hyper-luminous IR galax- 
ies). Overall, the large scatter in this diagram, and 
the strong dependence on the UV-slope of the un- 
known extinction curve serve as ca veats in apply- 
ing st andard UV dust-corrections. iBoquien et al. 



(2009) caution that applying the wrong extinction 
law to the UV would results in over- or under- es- 
timating the true SFR by factors of up to ~ 7. We 
return to this point in the following section. 

4.8. "IR Excess" and Obscured AGN Lu- 
minosity 

As discussed in § I4.7[ the BzK method has been 
used to study 24 /im galaxies in the GOODS field 
(|Daddi et al l l2007h . This paper found that as 
much as (20-30)% of 1.4<z<2.5 BzK selected 
galaxies with Ky e ga ^ 22 have infrared excess, in 

the Sense that SFRjr + SFRtjv, uncorrected 

>3.16x 

SFRtjv corrected- 

Using the X-ray stacking method 
on these 24 /am galaxies with "IR excess" , this 



study found that these galaxies contain highly ob- 
scured, Compton thick AGNs, and the predicted 
unobscured X-ray luminosity L 2 -8 keV ~ (1 -4) x 10 43 ergs/s. 
All our z ~ 2 ULIRGs are BzK galaxies satisfying 
the IR excess criteria (see Figure 120]) . 

In principle, if we know the correct relationship 
to estimate Ljr from rest-frame Lg^m, the correct 
UV extinction, and if the AGN contribution is neg- 
ligible, we should get SFR uv ^ corrected = SFR IR + 
SFRuv, uncorrected- The fact that all our z~2 
ULIRGs have "IR excess" could be due to either 
(a) non-negligible AGN contributing to total lu- 
minosities , or (b) wrong UV extinction correc- 



tion, o r (c) incorrect estimate of Lir. iDaddi et al 



( 20071 ) has argued that this "IR excess" is primar- 
ily due to Compton-thick AGN heating of mid-IR 
luminosities, causing the over-estimation of L/^, 
thus "excess" of SFR IR + SFR uv , uncorrected . Our 
mid-IR spectra suggest that most z ~ 2 ULIRGs 
(92%) have strong PAH emission, and AGN con- 
tribution to their Lir is small. However, the rest- 
frame Lsam estimated from observed broad band 
24 /am photometry will be over-estimated because 
of broad, strong PA H emission among s tarbursts. 
As already note by iMurphv et al. (|2QQ9h , the re- 
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Fig. 21. — Ratio of far-IR to UV flux at 1600 A compared to the UV spectral slope B for our galaxies an d some published 
comparison galaxies. The solid line shows the linea r fit to the A ifinn relationship by Meurer et al. ( 1999h. The dashed line 
show the relati on for normal star-f orming galaxies (|Cortese et al.ll2006l V The samples from Meure r et al . (1999); Reddv et al. 
(2006), and Chapma n et al.l (|2005h are marked with green triangles, cyan dots, and yellow squares, respectively. Our z ~ 1 
and z ~ 2 sample are marked with red and blue diamonds, respectively. The typical error of the points is shown in the bottom 
right corner of the plot. Most of our z ~ 2 ULIRG lie above the curve, in the region occupied by SMGs and HLIRGs. 



lationship used by iDaddi et al.l (|2QQ7h to com- 
pute Lir from Ls^m over-estimates Lir. For our 
ULIRGs, the median ratio between SFR estimated 
from Ls^m and is 3.1. Lastly, another 

big factor causing this "IR excess" is the under- 
estimated UV dust extinction correction. This is 
further supported by Figure EH where the z ~ 2 
ULIRGs are systematically offset above the local 
relationship defined by low luminosity, low extinc- 
tion star- forming galaxies. It is not surprising that 



UV dust extinction of heavily obscured ULIRGs is 
only a lower limit because, for a mixed distribution 
of stars and dust, increasing amounts of dust will 
further absorb stellar emission, but the reddening 
of the integrated colors of a galaxy will saturate as 
the most obscured spectral regions have less and 



less w eight in the integrated SED (jWuyts et al 



2009bD . 



Our IRS spectra illustrate that the "IR excess" 
feature in BzK galaxies is really due to strong 
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Fig. 22. — Fraction of the mid-infrared emission due to the AGN component and its X-ray luminosity (top and bottom 
panels, respectively) for the LIRGs and ULIRGs of our sample whi ch are not d ominated by AGN emission. The grey band 
marks the expected X-ray luminosity of obscured AGN according to Dad di et al . ( 2003). 



broad PAH emission, particularly 7.7 /im feature, 
produced by recent starbursts. This excess is not 
due to dust being heate d by obscured, Com pton 
thick AGNs as stated in iDaddi et"aH (12007). Be- 
low, we will use the IRS spectra to qualitatively 
estimate obscured AGN luminosities. 

Many mid-IR spectroscopic studies have shown 
that galaxies with obscured AGN tend to have 
excess mid-IR continuum at 4-6 /im, which is 



thought to be produced by small, h ot dust grains 



heate d by UV photons from AGN (jLaurent et al 
2000f ). There is a narrow correlation between the 
AGN luminosity at 5.8 /im, Ls.g/xm? and its intrin- 
sic X-ray luminosity L2-iokeV for a local samples 
of Seyferts and QSOs: 

logL 5 . 8/im = 1.209 logL 2 -io/eeV - 8.667, (3) 



with luminosities expressed in erg/s (jLutz et al 
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20041 : iBauer et al.ll2009h . 

To estimate the continuum due to the AGN 
emission at 5.8/im, we fit to each spectrum a com- 
bination of a linear component pl us the scaled 
averag e starburst template from iBrandl et al. 
(|2006l ) considering for the fit the 5.5-6.85/im in- 
terval. This direct approa ch provides r esult s con- 
sistent to the approa ch of iLutz et al] (|2004j ) and 
IValiante et all {2009) who used M82 as starburst 
template and adjusted for differences between M82 
and an average starburst g alaxy. To com pare our 
results to the work of Daddi et al. 1 (l2f)f)7h . we fol- 
low his criteria excluding galaxies detected in the 
2-8keV band in the 1 Ms Chandra survey and 
those without a clear decreasing in their SED af- 
ter 1.6 /im. Therefore, we do not consider the two 
ULIRGs whose mid-IR spectra show very little 
PAH emission, have a power-law SED and are 
clearly dominated by nuclear emission. We also 
exclude U4958 since its optical spectrum shows 
the presence of an AGN and the mid-IR spectra is 
almost featureless. Since several of the LIRGs of 
our sample are also "IR excess" galaxies, we apply 
the same analysis to this sample excluding three 
objects (L3945, L5511, L6211) which are not in 
the correct redshift range or have little or no PAH 
emission in their spectra. 

To estimate the AGN contribution to the to- 

r 



tal mi d-IR emission, we proceed as iMurphv et al. 
(|2Q09h by scaling the SED of Mrk231 to the 5.8 /mi 
AGN continuum and fitting the residual SEP wit h 
the average starburst from IBrandl et al.l (j2QQ6h . 
As shown in Figure [22j the emission from most of 
the ULIRGs and almost all LIRGs is completely 
dominated by star for mation. This c onfirm s and 
extends the analysis of lMurphv et al.l (|2009f ) , since 
our study is done on a bigger, fainter, and more 
homogeneous sample. Moreover, our analysis al- 
lows us to directly estimate the X-ray luminosity 
of the AGN harbored in the galaxies of our sam- 
ple. In the bottom panel of the same Figure [22j 
we report the estimated X-ray luminosities of the 
AGN component of the LIRGs and ULIRGs of 
our sample compared to the predictions b a sed o n 
stacked X-ray analysis from iDaddi et al.l ([20071 ). 



Also in this case, the majority of ULIRGs have 
AGN with X-ray luminosities lower than these val- 
ues. A straight average of the values found would 
give a value of L2-iokeV — (1.8±0.6) x 10 43 erg/s 
which would be compatible with the analy- 



sis of IDaddi et al.l ([20071 ). A robust average, 
using a biweight estimator, gives a value of 
£2-10 kev = (0.1±0.6) x 10 43 erg/s. This means 
that a few galaxies with bright hidden AGN dom- 
inate the stacking of the faint infrared sources in 
the aforementioned analysis, highlighting the in- 
trinsic danger of any stacking analysis. We can 
conclude that the infrared emission from the ma- 
jority of our ULIRGs (21/24, i.e. 88% of our 
sample) is dominated by star formation and the 
analysis of our sample of ULIRGs sets a low upper 
limit to the presence of Compton thick AGN in 
luminous infrare d galaxies at z ~ 2, with respect 
to the analysis of IDaddi et al.l (|2007h . 

In the case of LIRGs, the straight and robust 
average are (0.2±0.1) x 10 43 erg/s and (O.OiO.l) x 10 4 
respectively. We conclude that for the z~l sam- 
ple, the majority of the galaxies do not harbor 
hidden AGN and their emission is dominated by 
star formation. 

4.9. Morphologies 

LIRGs and ULIRGs shine at infrared luminosi- 
ties of 10 44 — 10 46 ergs/s. It is not yet clear what 
physical processes light up a galaxy in the infrared 
and if major mergers are the dominant mecha- 
nism at z ~ 1 — 2, as it is the case for local 
ULIR Gs ([Sanders fc Mirabellll996l : IVeilleux et al 



20021 ). The simplest way to address this question 
is to examine the morphologies of our sources us- 
ing the high spatial resolution images from HST. 

Figures [23] and [24] show image cutouts for 
ULIRGs and LIRGs, respectively, in the HST 
F435W (B), F606W (V) , F775W (i), F850LP 
(z) bands and the ground based J and K images. 
Each cutout has a size of 10 x 10 . These two 
figures illustrate the steep changes of the spectral 
energy distribution from B-K, and clearly indi- 
cate that many of our sources have extremely red 
optical-to-NIR colors. To better show colors and 
morphologies, we present three-color images of the 
galaxies which have been observed with ACS in 
the F435W, F606W and F850LP bands. The an- 
gular size of the images is 3 and 4 for ULIRGs 
and LIRGs, respectively, which corresponds to 
25.2 kpc and 31.2 kpc at the median redshift of 
1.9 and 0.9, respectively. To obtain the three- 
color images, we filtered each ACS image with a 
wavelet transform to get rid of the background 
noise and to enhance the faint and extended mor- 
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Fig. 23. — The image cutouts for each ULIRG are in the filters F435W, F606W, F775W, F850LP from HST ACS images of 
the GOODS treasury survey. J and Ks images are from the public EIS survey done with ISAAC on the VLT telescopes. The 
object ID is labeled on the F435W image. Each image stamp is 10 x 10 . 



phological features. The technique used to filter 
the signal is a 2-dimension generalization of th e 
method described in detail in lFadda et al.l ([ 19981 ). 
To each band we assigned the same weight. 



Many extended sources, such as ULIRG4367, 
4451, 5652 and 5829, are also very red, suggest- 
ing that dust obscuration could be in play. Only 
about 4 sources (428,4950,4958,5877) have very 
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Fig. 24. — Similar to Figure [231 this plot shows the image cutouts for 21 LIRGs at z ~ 1. We used all of the available 
observations in HST/ACS F606W and F850LP filters for LIRGs within the GEMS survey, and F435W, F606W, F775W and 
F850LP for LIRGs within the GOODS region. For objects in the GOODS region, J and Ks images are from the EIS survey. 
Otherwise, they are from the MUSYC survey (|Gawiser et al.ll2C)06h . 



compact or nearly unresolved morphologies. It is 
interesting to note that two of these compact mor- 
phology sources (ULIRG4958 and 5877) show very 
blue compact cores in Figure [25j and are AGN 
dominated systems with weak or no PAH emis- 
sion shown in Figure [7] This correspondence be- 
tween the rest-frame optical compact morpholo- 



gies and the mid-IR AGN dominant spectra has 
also been found among a large sample of 24 /im 
selected ULIRGs in the Spitzer First Look Survey 
(jZamoiski et al.ll2009l foasvra et al.ll2008l ). 

To better quantify morphological differences 
between our sample of galaxies and other galax- 
ies in the field which are faint infrared emitters, 
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Fig. 25. — Three color images for 24 ULIRGs and 10 LIRGs of our sample with F435W, F606W, and F850LP ACS coverage. 
The field of view is 3" x 3 and 4 x 4 for ULIRGs and LIRGs, respectively. The composed image was obtained by giving 
the same weight to each band. 



we computed the Gini coefficient and M20 for the 
galaxies in our sample and control samples of faint 
infrared galaxies in the same redshift ranges. The 



Gini coefficient, G, is a statistic to measure the dis- 
tribution of flux within the galaxy and M20 is the 
second-order moment of the brightest 20% of the 
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Fig. 26. — Gini coefficient and M20 statistic for the sample of LIRGs (red filled dots) at z~l, ULIRGs (blue filled dots) 
at z~2 and faint infrared sources in the redshift range z~l and z~2 (red and blue empty dots, respectively). The triangles 
correspond to the sources in our sample whose mid-IR emission is AGN-dominated. The grey shaded region highlights the 
Gini-M2o sequenc e which, locally, is populated by "normal" galaxies while local ULIRGs lie above the sequence. The three 
regions defined in lLotz et al. I (2008) are marked with solid lines. The statistics have been computed from ACS images in the 
bands F606W and F850LP for galaxies at z~l and z~2, respectively, to be approximately at the rest-frame wavelength of 
3100 A. 



galaxy's flux. These tw o nonparametric statistics 
have b een introduced bv lLotz. Primack fc Madaul 
( 2004 ) to quantify galaxy morphology. Normal 
galaxies form a sequence in the G— M20 plot and, 
in the local Universe, ULIRGs lie above this se- 
quence. To compute th ese statistics, we followed 
the m ethod described in lLotz. Primack fc Madau 
(|2004j ) using our wavelet-filtered images to reduce 
the role of noise in the computation. We computed 
G and M20 using the F850LP and F606W images 



for ULIRGs at z ~ 2 and LIRGs at z ~ 1, re- 
spectively, to be roughly at the same 3100 A rest- 
frame wavelength. We note that also the pixel 
size of the galaxies at z ~ 1 and z ~ 2 correspond 
to a similar physical scale (0.24 and 0.25 kpc, re- 
spectively, with our adopted cosmology), so that 
we can compare the values computed for the two 
redshift ranges. To compute the errors on the co- 
efficients, we proceeded using Monte Carlo simu- 
lations. It is, in fact, impossible to have a direct 
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Fig. 27. — The plot shows the HST/NICMOS H-band 3" x3" images centered on 10 ULIRGs and one LIRG of our sample. 
Labels correspond to the object IDs. 



error propagation since the computation of the co- 
efficients involves filtering of the image, selection 
of regions, sorting of values and search of source 
center. Therefore, for each source, we created 50 
synthetic images by adding a random realization 
of a noise image (assuming a Gaussian distribution 
with dispersion equal to the rms of the image) to 
the wavelet-filtered image. After adding the orig- 
inal zero, we rerun the code to measure the two 
coefficients on the synthetic images and estimated 
the errors assuming that they are equal to the dis- 
persion of the values computed from the synthetic 
images. The measured errors are 0.014±0.011 and 
0.04 ± 0.04 for the G and M20 coefficients, respec- 
tively. 

In Figure [26j we report the G and M20 values 
for the galaxies of our sample for which the signal- 
to-noise was sufficient to compute them and values 
for a control sample of galaxies with faint infrared 
emission (less than 20/iJy at 24/im) in the redshift 
ranges 0.8-1 and 1.8-2 where we find most of the 
galaxies of our sample. The typical l-cr error is 
reported in the lower right corner. A grey shaded 
region highlights the G-M20 sequence which, lo- 
cally, is typically populated by "normal" galaxies 
while ULIRGs are found above that seq uence. The 



three regions of the plot defined by lLotz et al. 
(2008) populated by merger, early-type and late- 
type galaxies are marke d with solid l ines. The G- 
M 20 sequence found by lLotz et al. I |2008) in the 
Groth strip at 0.2 < z < 1.2 in the rest-frame blue 
band is well reproduced by our study, although 
at a bluer rest-frame wavelength. We immedi- 



ately notice that the separation between infrared- 
active and infrared-quiet galaxies is not so clear as 
for normal and IR- luminous local galaxies. Most 
of the galaxies of our sample follow the same G- 
M20 sequ ence as the majority of the infrared- quiet 



galaxies. lLotz et al. I (|2008[ ) found that the per- 
centage of 24/im sources at 0.4 < z < 1.0 with 
Lir > 10 n L© are disk galaxies and only ~15% 
are classified as major merger candidates. At the 
3100 A rest-frame, our study shows that z ~ 1 
LIRGs have a higher merger rate (37±^%) and 
that only 21+^% of the ULIRGs at z ~ 2 can be 
classified as mergers. We can therefore conclude 
that the morphologies of most of the IR-luminous 
galaxies at z ~ 1 and z ~ 2 do not differ in most 
of the cases from those of "normal" infrared-quiet 
galaxies. Repeating the analysis for the LIRGs 
of our sample in the blue rest-frame, i.e. using 
the F850LP band ACS images, only two galaxies 
populate the "merger" region for a percentage of 
10^7 4 %, w hich is perfectly com patible with the re- 



suits from lLotz et al. I (2008). This result is not 
surprising since the rest-frame 3100 A light is more 
sensitive to massive, young stars so that late- type 
galaxies appear more patchy and have no promi- 
nent opti cal bulge with re s pect to optical images 

In the case of 



(see, e.g. iKuchinski et al. 



s pect to 
([19991)). 



ULIRGs, we have HST NICMOS images in the 
H band, approximately corresponding to the blue 
rest-frame at z ~ 1.9, for a subset of 10 ULIRGs 
(see Figure [27j) . The Gini coefficients have a typi- 
cal difference of 0.02 with respect those computed 
in the ACS images. All these sources are classified 
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as "normal" galaxies in the G-M20 plot. Unfortu- 
nately, we do not have NICMOS coverage for the 
three ULIRGs which fall in the merger region of 
the G-M20 plot. 

5. Summary 

The study presented in this paper aims to char- 
acterize the mid-IR spectral properties of very 
faint 24 /im sources with SW/zm ~ 0.15-0.45 mJy, 
specifically targeting LIRGs at z ~ 1 and ULIRGs 
at z ~ 2. This sample contains the faintest sources 
observed with IRS, reaching the limit of the in- 
strument capability. These 24 /im faint LIRGs and 
ULIRGs are infrared luminous, sampling the typi- 
cal luminosities of infrared luminosity functions at 
z ~ 1 & 2. Furthermore, they are similar to normal 
galaxies selected by optical/near-IR surveys in the 
sense that their stellar masses are < M*, the turn- 
over of the mass functions (~2x 1O 11 M ). The 
mid-IR spectra, in combination with X-ray, IR and 
optical observations, have allowed us to determine 
physical properties of high-z LIRGs and ULIRGs, 
putting them in the context of larger galaxy pop- 
ulations selected at optical/near-IR wavelength. 
The major results of this work can be summarized 
as follows: 

• The vast majority of our sample (45 out 
of 48) has IRS spectra with clear PAH emission 
and/or silicate absorption, allowing secure redshift 
measurements. In particular, a third of the z ~ 2 
ULIRGs are DRGs (galaxies with (J - K) > 2.3) 
for which redshifts are very difficult to determine 
from optical or near-IR spectroscopy. Our finding 
of DRGs being dusty star forming galaxies, rather 
than passively evolving old stellar population is 
consistent with previou s studies bv (jLabbe et al.l 
l2QQ5t IWuvts et al.ll2QQ7t IPapovich et al.ll2006f ).~ 

• The fraction of AGN harbored by the sources 
in our sample is much smaller than that of other 
samples of more luminous distant ULIRGs in lit- 
erature. From the analysis of our mid-IR sources, 
we estimate that their total IR emission is mainly 
powered by star formation. AGN emission dom- 
inates the energy output of LIRGs and ULIRGs 
only for 5% and 12% of them, respectively. 

• A quantitative analysis of the morphology 
based on ACS images reveals that most of the IR- 
luminous galaxies of our sample do not differ sig- 
nificantly from the majority of the IR-quiet galax- 



ies at the same redshift range. At the 3100 A rest- 
frame, 37±?|% of the z ~ 1 LIRGs and 21±^% 
of the z ~ 2 ULIRGs have Gini coefficients higher 
than those of normal IR-quiet galaxies. In the blue 
rest-frame, 10ty 4 % of the z ~ 1 LIRGs fall in the 
merger region of the Gini-M2o plot which is com - 
patible with the 15% found bv Lotz et al~l (|2008h 
among infrared galaxies with Lir > 10 n L© at 
0.4 < z < 1.0. 

• We found a maximum value of ~ 1 jam for 
2/xm for all infrared galaxies, independent of 

redshift and luminosity. This ceiling in 6.2 /im 
PAH strength is due to the fact that, for pure 
starbursts, the strength of the PAH line and of 
the underlying continuum are both proportional 
to the luminosity due to star-forming activity. At 
^ 10 n L©, there is a clear anti-correlation 
between the rest-frame 24 /im luminosity and 
6.2 /im PAH equivalent width. The 6.2 /im PAH 
equivalent width of the z>l sources is systemati- 
cally greater than that of the local luminous galax- 
ies in the same luminosity range. This is likely due 
to an AGN contribution in high-redshift galaxies 
lower than that of local IR luminous galaxies. 

• All our z ~ 2 ULIRGs ar e BzK galaxi e s with 



infrared excess, according to iDaddi et al.l <J_2J 



definitions. Stacking all the sources whose infrared 
flux is not dominated by AGN emission, we pro- 
vide quantitative constraints on the dust obscured 
AGN emission at the rest-frame 6 /im. By sub- 
tracting from the spectra around the 6.2 /im PAH 
feature the scaled averaged starburst spectrum, we 
measured the average AGN powered continuum at 
6 /im, estimating an average intrinsic X-ray AGN 
luminosity of L2-iokeV = (0.1 ±0.6) x 10 43 erg/s, 
a value substantially lower tha n that inferr e d from 
the X-ray stacking analysis in IDaddi et al. I (j2007h . 
A similar analysis on the LIRG sample suggests 
that their emission is mostly powered by star for- 
mation. 

In conclusion, our mid-IR spectra of LIRGs and 
ULIRGs at z ~ 1 and 2 reveal that the majority of 
these galaxies have strong PAH emission and are 
starburst dominated. The Compton thick AGN 
contribution to bolometric luminosity for this sam- 
ple selected with 624/1™ °° 0.15-0.5 mJy is small. 
The so-called "IR excess" for these faint 24 /im se- 
lected ULIRGs at z ~ 2 is mainly due to the strong 
PAH contribution to the rest-frame mid-IR lumi- 
nosity which leads to over-estimations of Lir and 
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to under-estimation of the UV dust extinction. 
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Fig. 28. — Analysis of an IRS observation through the GUI of the IrsLow code. In this case, the target is a double source. At 
least other three serendipitous sources are visible. The two images show the 2nd and 1st orders after coadding all the frames 
(including the parallel observations). In the right upper panel, the profiles of the averaged rows of the two orders are displayed. 
The shaded regions have been masked to compute better background and rms images. Spectra of the two sources between 
the green lines on the images have been extracted and are shown in the lower left panels. Finally, the PSF spectral fit at the 
wavelength corresponding to the vertical line is shown in the right-bottom window. Each point corresponds to a pixel in the 
original frames. Circled points are masked for the fit. The lower boxes contains the limits of the image and coordinates of the 
extracted spectrum. 

A. IrsLow, a software package for low-resolution IRS spectra 

For the analysis of the spectra presented in this paper, we developed a new software package called 
"IrsLow" . This software allows one to exploit in an optimal way the redundancy of the data taken with the 
IRS spectrograph and to extract spectra in a very flexible way. The software and a tutorial are available 
onlinqj. The code requires as input the basic calibrated data (BCD) produced by the Spitzer pipeline 
and uses the full redundancy of the data to reject bad pixels and compute the residual background of the 
observations. Once the bad pixels are masked and the background subtracted from each individual image, 
it coadds all the images of each order at different positions. The parallel data of non-primary orders are 
coadded to identify and extract spectra of serendipitous sources. The spectra are extracted using PSF fitting 



web.ipac.caltech.edu/staff/fadda/IrsLow 
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Fig. 29. — Spatial PSF profiles (left) and FWHM as a function of the wavelength (right) in the different orders computed 
with calibration stars. Amplitude and FWHM of the PSF profiles are normalized to one. These PSFs are used for the optimal 
extraction with IrsLow. 



and the code allows the simultaneous extraction of two spatially close spectra. An example output from this 
code is shown in Figure [281 

The optimal background subtraction is done iteratively with manually masking off all spectra. Uncertainty 
and background images are obtained fro m the stacked BCD after the rejection of pixels affected by cosmic 
rays or sources. A biweight estimator (|Tukevl[l958f ) is used to compute robust estimates of average and 
dispersion also with a low number (~ 10) of frames. After extracting the spectrum, each PSF fit for a single 
wavelength position can be visualized to delete other possible outliers and therefore improve the fit. Pixels 
with highly varying responsivity (i.e. rogue pixels in the IRS jargon) are identified by comparing their RMS 
values to those of the neighboring pixels. The final image is obtained by fitting a B-spline at each wavelength 
to all the flux values from different pixels at their relative spatial positions. The fit is repeated three times, 
discarding 7-cr and 4-<j outliers in the second and third iteration, respectively. We derive PSF profiles using 
the IRS calibration stars. The average profiles are obtained by normalizing the flux of the spectra and scaling 
the FWHM as a function of the wavelength. The scaling of the FWHM is shown in the same Figure [29j 
Finally, the extracted spectra are saved in ASCII files with names equal to their spatial coordinates. The flux 
and wavelength calibration files are the same as what is included in the IRS software SPICE. The calibration 
version is automatically selected by the software according to the information stored in the BCD headers. 
In fact, the calibration version is linked to the version of the pipeline used to obtain the BCD files. 

One unique feature of our IRS reduction software is how it handles confused spectra from two close-by 
sources. As shown in Figure [281 t ne software allows the optimal extraction of two close sources at the same 
time by using our derived PSF profiles, shown in Figure [29j This feature has been particularly useful in 
the case of the source U16526. The observed IRS spectrum for U16526 is, in fact, a blend of two spectra 
at z = 0.37 and z— 1.77 (photometric redshifts from the FIREWORKS catalog), spatially separated by only 
6 . This corresponds to the spatial resolution of the IRS spectrograph in the wavelength range used in the 
observation. When displaying this spectrum with our software, one can see only that the major features of 
the two sources are slighlty spatially displaced although the continuum of the two sources is blended. In our 
first extraction, we erroneously considered one single source and interpreted some strong emission features 
as H2 lines. Then, after examining the 24/im image of the source, we realized that the spectrum was a blend 
of spectra from two different sources at different redshifts. The unique deblending capability of our software 
allowed us to separate the spectra of these two confused sources. In Figure [30j the left panel shows the total 



39 




12 15 18 21 24 



A obs/( 1+Z ) t/ 2 " 1 ] 

Fig. 30. — Spectral extraction of ULIRG 16526 considering one object (left) and fitting two objects (right). If only a single 
spectrum is extracted, spectral features from a nearby low-redshift object (at z = 0.37) could be wrongly interpreted as strong 
H2 and 6.2 |Um PAH features. 

spectrum, which is the sum of two spectra, one at z = 0.37 and one at z=1.75, illustrated on the top right 
and bottom right, respectively. At this point it was clear that the H2 lines in the blended spectrum were in 
fact only strong emission features of the z = 0.37 source. 
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Table 2 

The Sample Targets and IRS Observation Log 
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a SL and LL stand for Short-Low and Long-Low modules. SL 1st order, LL 2nd order, and LL 1st order 
cover the ranges 7 — 14/ira, 14 — 21/ira, and 21 — 35/ira, respectively. 
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Table 3 

The Sample Target Fluxes and Redshifts 



ID 






um 




urn 




ZIRS 




* 


log(L 2 4/^m) 


M* 


EW6.2fj,m 




[A*Jy] 


SN 




SN 




SN 


[raJy] 


SN 








[L©] 


[1O 1O M ] 


[H 


U428 


15.9 


68 


65 


24 


287 


38 


2.43 


3 


1.783 


(1.664) 


a 


11.47 


4.09 


1 on_Ln OA 

l.o9±U.o4 


U4367 


22.2 


79 


101 


12 


169 


20 


<1.55 




1.624 


(1.762) 


a 


11.23 


13.10 


0.54±0.z3 


U4451 


11.1 


38 






225 


28 


<1.53 




1.875 


(1.684) 


a 


11.40 


3.47 


0.83±0.50 


J J A a on 

U4499 


16.3 


77 


< 35 




182 


23 


<1.45 




1.956 


(1.909) 


a 


11.58 


11.90 


>0.84 


U4631 


15.4 


60 


< 21 




269 


33 


<1.40 




1.841 


1.896 


1 


11.39 


9.61 


>0.92 


U4639 


13.5 


57 


< 14 




216 


27 


<1.45 




2.112 


2.130 


5 


11.58 


3.57 


>2.61 


U4642 


12.3 


60 


65 


7 


268 


32 


2.14 


3 


1.898 


(1.748) 


a 


11.82 


2.50 


1.33±0.41 


U4812 


15.9 


42 


< 53 




317 


41 


3.75 


13 


1.930 


1.910 


1 


12.00 


6.63 


>1.04 


U4950 


163.3 


598 


453 


54 


547 


67 


<1.68 




2.312 


2.291 


4 


11.79 


10.41 


r> i o I r> r\r' 

0.1z±0.0o 


U4958 


16.4 


62 


109 


11 


262 


33 


<1.67 




2.118 


2.145 


5 


11.51 


8.50 




U5050 


24.6 


97 


51 


7 


194 


22 


<1.50 




1.938 


(1.726) 


a 


11.40 


18.93 


O.o0±0.34 


U5059 


18.1 


64 


< 37 




253 


31 


<1.53 




1.769 


(1.543) 


a 


11.49 


6.59 


>1.61 


U5150 


13.1 


35 


65 


6 


286 


37 


2.98 


12 


1.898 


(1.738) 


a 


11.92 


2.17 


U.77±0.47 


U5152 


13.1 


43 


66 


7 


305 


38 


<1.77 




1.794 


(1.888) 


a 


11.44 


2.75 


0.50±0.34 


U5153 


19.2 


68 


< 34 




183 


21 


<1.71 




2.442 


(2.030) 


b 


11.73 




>0.41 


U5632 


19.2 


86 


83 


35 


421 


52 


<1.38 




2.016 


1.998 


1 


11.75 


6.79 


0.99±0. 22 


U5652 


30.3 


153 


154 


17 


316 


44 


3.24 


8 


1.618 


1.616 


1 


11.68 


15.73 


l.zl±0.z3 


U5775 


9.6 


51 


< 31 




141 


16 


2.64 


11 


1.897 


(1.779) 


a 


11.61 


5.30 


>0.72 


U5795 


18.7 


63 


< 31 




264 


33 


<1.43 




1.703 


(1.524) 


a 


11.36 


5.63 


>2.16 


U5801 


16.2 


57 


< 35 




216 


24 


<1.45 




1.841 


(1.642) 


a 


11.25 


13.68 


>0.83 


U5805 


8.0 


25 


68 


11 


170 


19 


<1.45 




2.073 


(2.093) 


a 


11.74 


1.16 


0.41±0.18 


U5829 


15.4 


69 


57 


12 


208 


26 


<1.40 




1.742 


(1.597) 


a 


11.32 


9.21 


0.98±0.31 


U5877 


31.7 


120 


179 


23 


360 


48 


3.88 


19 


1.886 


(1.708) 


a 


11.77 


3.54 


0.17±0.10 


T T 1 C C 0£? 

U loozo 


20.6 


66 


69 


18 


306 


45 


5.17 


13 


1.749 


(1.718) 


a 


12.00 


4.05 


0.53±U.35 


L3832 


23.5 


61 






207 


28 


<1.90 




0.767 


(0.763) 


b 


10.53 


4.39 


0.81±0.4o 


L3945 


19.0 


18 






367 


39 


3.81 


11 




(0.950) 


b 


11.08 


0.37 


0.20±0.05 


L4020 


20.8 


30 






332 


38 


4.22 


21 


0.826 


(0.779) 


b 


10.90 


3.10 


1.13±0.13 


L4177 


20.5 


70 


204 


22 


208 


27 


<1.78 




0.842 


0.840 


3 


10.62 


3.71 


0.91±0.2o 


L4419 


46.8 


159 






259 


38 


4.06 


13 


0.974 


(0.995) 


b 


11.08 


11.23 


r\ An I n no 

0.4o±0.z3 


L4900 


26.8 


95 


267 


32 


213 


31 


<1.82 




1.047 


1.045 


5 


10.91 


6.75 


1.11±0.05 




22.6 


106 


9ftQ 
ZObJ 


1 1 n 
± ±u 


99Q 

zzy 


Zo 


<1.42 




i.uoy 


l.UoO 


i 
i 




^ 09 


i Ozi-i-n ^9 


L5367 


24.3 


27 






252 


26 


<1.92 




0.974 


(0.868) 


b 


10.86 


5.79 


1.06±0.07 


L5420 


36.7 


135 


450 


48 


381 


46 


2.76 


5 


1.068 


1.068 


3 


11.15 


4.42 


0.91±0.20 


L5511 


38.6 


50 






401 


54 


2.22 


5 




( 2.08) 


c 


11.87 


25.90 


0.48±0.23 


L5630 


22.8 


104 


277 


149 


286 


36 


2.06 


3 


0.997 


0.996 


3 


10.93 


2.27 


1.24±0.16 


L5659 


28.6 


129 


285 


25 


231 


30 


2.41 


9 


1.044 


1.038 


1 


11.00 


4.57 


1.13±0.53 


L5876 


29.2 


103 


310 


41 


330 


43 


2.54 


15 


0.971 


0.969 


1 


10.98 


7.16 


1.11±0.04 


L6211 


31.5 


84 






457 


62 


<1.88 




1.843 


(1.840) 


c 


11.70 




1.56±0.22 


L6221 


19.1 


14 






239 


12 


<2.87 




1.012 


(0.989) 


b 


11.02 




1.04±0.26 


L6473 


32.3 


44 






484 


63 


2.78 


8 


0.816 


0.811 


1 


10.83 




1.10±0.28 


L6658 


45.2 


60 






467 


59 


<1.87 




0.969 


(0.963) 


b 


10.94 


12.54 


1.33±0.24 


L7079 


23.3 


20 






289 


27 


3.21 


8 


0.955 


(0.884) 


b 


11.00 




1.01±0.01 


L13667 


33.8 


30 






496 


64 


<2.65 




0.936 


(0.924) 


b 


11.00 




1.57±0.01 


L13958 


21.8 


83 


292 


39 


376 


47 


3.31 


13 


0.891 


0.896 


2 


10.96 


2.10 


1.99±0.24 


L14143 


34.8 


22 






427 


35 


4.71 


8 


1.043 


(0.965) 


b 


11.29 




0.89±0.08 


L14228 


23.8 


23 






313 


35 


<2.12 




0.953 


(0.961) 


b 


10.89 


3.02 


1.55±0.21 


L15906 


18.4 


63 


193 


23 


262 


41 


3.23 


8 


0.976 


(1.045) 


a 


11.02 


0.44 


1.13±0.17 


L16005 


24.2 


32 






253 


42 


3.29 


9 


0.865 


(0.984) 


b 


10.86 


1.59 


1.60±0.06 



* Sp ectro scopic redshifts from : (1) IVanzella et al.l (|2008l . l2006h . f2)lMignoli et al.l (120051 ). f3)lLe Fevre et al.l (120041) . (4^ ISzokolv et al.l 
(2004), f5) lPopesso et al.l (120081) . Photometric redshifts from: (a) IWuvts et alJ (l2008T )77b) ICaputi et al.l (120041 ). (c) this paper. 
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